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ABSTRACT 

In the z = 1.94 Civ absorption line system in the spectrum of quasar Q1222+228 
(z em = 2.04), we find two clouds which have contrasting physical conditions, although 
they are only at a 17 km s~ 1 velocity separation. In the first cloud Sin, Siiv, and 
Cn are detected, and Alii and Aim column density limits in conjunction with 
photoionization models allow us to infer that this cloud has a large Si abundance and a 
small Al abundance relative to a solar abundance pattern. This pattern resembles that 
of Galactic metal-poor halo stars, which must have formed from such high redshift 
gas. The second cloud, in contrast, has detected Alii and Aim (also Siiv and Cn), 
but no detected Sin. We demonstrate, using photoionization models, that Al/Si must 
be greater than (Al/Si)© in this unusual cloud. Such a ratio is not found in absorption 
profiles looking through Milky Way gas. It cannot be explained by dust depletion 
since Al depletes more severely than Si. Comparing to other Al-rich environments, we 
speculate about the processes and conditions that could give rise to this abundance 
pattern. 



Subject headings: quasars: absorption lines — galaxies: structure — galaxies: evolution 



Introduction 



One of the motives for studying quasar (QSO) absorption line systems is to document the 
photoionization structure, chemical composition, and kinematics of the ISM and halos of high 
redshift galaxies at a level of detail on par with studies of the Milky Way ISM and Halo. QSO 
lines of sight that pass through high redshift galaxies sample multiple gaseous structures having a 
variety of physical conditions such as found in Galactic Hi and Hn regions, supershells, infalling 
Halo gas, and material being processed at the Galaxy/Halo interface. Recognizing the abundance 



1 Based in part on observations obtained at the W. M. 
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patterns, photoionization conditions, and kinematics associated with these various types of 
structures is a prerequisite for a detailed understanding of the evolution of galactic gas. 

Based upon low resolution spectra, early QSO absorption line efforts have been limited to 
simple curve of growth analyses of equivalent widths, providing only a weighted-average of the 
absorbing gas properties in each system. A great deal has been learned from these studies and 
a global (statistically based) picture of chemical and ionization evolution has been suggested 
[see Steidel (1993) and references therein]. The next logical step is to examine variations in the 
chemical and ionization conditions within single galaxies and to incorporate kinematic information. 
Using high resolution UV spectra, several researchers studied the Milky Way Disk and Halo at a 
detailed level ( Welty et al. 1997 ; Savage fe Sembach 199l| ; | now et al. 1996 { [Spitzer fe Fitzpatrick| 
1995| ; |Fitzpatrick fc Spitzer 1994] ; |Spitzer fe Fitzpatrick 1993| ). Some recent efforts have focused on 



the cloud-by-cloud conditions in high redshift galaxies (cf. [Tripp, Lu, fc Savage 1997 ; Petit jean 



et al. 1994 ). It is hoped that, ultimately, a statistical picture of the processes that give rise to the 
observed absorbing gas properties will improve our understanding of the present-epoch Milky Way 
and galactic evolution. In this Letter we study the cloud-to-cloud properties in a Civ system at 
z ~ 1.94 along the line of sight toward the quasar Q1222 + 228 in order to: 1) demonstrate a large 
variation in ionization and/or abundance conditions between two kinematically adjacent clouds in 
the same absorber; and 2) present an unusual cloud that has an Al/Si abundance ratio enhanced 
by a factor of several relative to the solar ratio. 



2. The Data 

Our spectrum of Q1222 + 228, which has resolution Aufwhm = 6.6 km s _1 with three pixels 



per resolution element, is a 3600 second exposure obtained with HIRES/Keck I ( Vogt et al. 1994 ) 
on 23 January 1995. The wavelength coverage is 3810.5-6304.9 A with breaks redward of 5100 A. 
The spectrum was obtained as part of a large study of intermediate redshift Mgn absorbers 



( Churchill 1997 ), and thus the captured wavelength coverage did not include certain desirable 
transitions of the z ~ 1.94 Civ system. The reduction and analysis of the data are described 
elsewhere (Churchill 1997, Churchill, Vogt, & Charlton 1998a, Schneider et al. 1993, Churchill 
et al. 1998b). In Figure |l|, we present the detected transitions (5a) aligned in line-of-sight velocity. 

We focus upon the Sill A1527, Alii A1671, and Aim A1863 transitions in two absorbing 
clouds centered at v ~ and ~ 17 km s . The Civ AA1548, 1550 profiles are saturated across 
this velocity interval, but the Siiv and Cll profiles clearly show two distinct clouds. In Cloud A, 
at v ~ km s , we have detected Sin A1527 and obtained upper limits on both Alii A1671 and 
Aim A1863. In Cloud B, at v ~ 17 km s , the converse is true; we have detected Alii A1671 and 
Aim A1863, and obtained an upper limit on Sin A1527. The Aim A1855 transition, the doublet 
counterpart of Aim A1863, was not covered. 



The cloud velocities, column densities, and Doppler b parameters were obtained from Voigt 
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profile (VP) decomposition of the spectra using the program MINFIT ( Churchill 1997| ), which 



minimizes x 2 between the model and the data. All detected transitions were fit simultaneously, 
with the exception of the Civ A1550 transition, which is blended with an unidentified feature. The 
VP results are listed in Table [l]. The 5<r equivalent width limits for the non-detected transitions 
were converted to column density upper limits using the curve of growth. Since these equivalent 
width limits correspond to the linear part of the curve of growth, the column density limits are 
not sensitive to the adopted b parameter (see Table [l]). We note that iV(Siii) also corresponds to 
the linear part of the curve of growth and therefore is well constrained despite the large fractional 
error in 6(Sill). 

Constraints are available for the total neutral hydrogen column density of the absorber, which 
can be interpreted as an upper limit on iV(Hi) of any single cloud in the system. In the FOS/HST 
spectrum of Impey et al. (1996), there is no indication of a break at the expected location of the 
Lyman limit. They report detections of Ly/3 and Ly7 that are members of complex blends which 
suffer from a pre-COSTAR point spread function. Fernandez-Soto et al. (1995) performed a 
profile fit to a saturated, multi-component Lya feature at ~ 3571 A in a spectrum with resolution 
~ 30 km s -1 . However, the number of components and their b parameters and column densities 
are too uncertain to provide a constraint more useful than that provided by the lack of a Lyman 
limit break. We have conservatively adopted an upper limit of log iV(Hi) ~ 16.5 cm -2 . 



3. Photoionization Modeling of the Clouds 

We inferred the clouds' physical conditions using CLOUDY ( [Ferland 1996j ). We have assumed 
a, z = 2 extragalactic UV ionizing background, given by Haardt & Madau (1996) (logF u = —20.35 
erg s _1 Hz -1 cm -2 ). The input parameters for CLOUDY are the neutral hydrogen column 
density, iV(Hi), the metallicity, Z, and the total number density of hydrogen, tih- CLOUDY 
assumes that photons are incident on one side of a plane-parallel slab of gas with constant n# 
and integrates the equations of radiative transfer along an optical path through the slab until 
it reaches a specified A r (Hi). For a fixed spectrum, there is a relationship between iih and the 
ionization parameter U (the number density of hydrogen ionizing photons over the number density 
of electrons). For simplicity, we adjusted the two ratios^] [Si/H] and [Al/H] keeping [C/H] = 0. An 
important assumption is that the modelled transitions arise in a single phase environment. The 
broad, saturated, Civ profile may arise in a physically distinct higher ionization phase from the 
Cn and other lower ionization species. We argue that the precise velocity alignment of the lower 
ionization species and the fact that their b parameters are consistent (see Table |]) justify the 
assumption. 

For Cloud A, we measured iV(Siii), iV(Siiv), iV(Cn) and upper limits on iV(Alii) and 



3 We use the notation [X/H] = log(X/H) - log(X/H) . 
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iV(AllIl). Without precise information on the neutral hydrogen distribution in the system, we 
modelled the cloud for a range of iV(Hi). For a given iV(Hi) we ran CLOUDY in its "optimized" 
mode, allowing Z, njj, and [Si/H] to vary until the cloud model was consistent with the observed 
iV(Siii), A(Siiv), and N(Cu). The ratio iV(Sin)/iV(Siiv) constrained both the ionization 
parameter and metallicity. For models with logiV(Hi) > 15.4 cm~ 2 , iV(Alii) was greater than the 
upper limit. We therefore set all other parameters equal to those given by the optimized solution 
and then solved for the model cloud by varying [Al/H] and optimizing on the measured upper 
limits on iV(Alll) and iV(Aliii). In the upper panel of Figure ^, we illustrate how the [Si/H] 
and [Al/H] ratios, U, and Z, must be scaled for the models to be consistent with the data. For 
log iV(Hl) > 15.5 cm -2 , a factor of two to three enhancement of [Si/H] and an underabundance 
of [Al/H] are obtained, nearly independent of metallicity. The lack of a Lyman limit break in 
the FOS/HST spectrum of Impey et al. (1996) precludes a large A(Hi). This cloud could be 
consistent with [Si/H] and [Al/H]~ (solar values) only if the metallicity is super-solar [for 
logA(Hi) ~ 15 cm -2 ]. We conclude that [Si/H] is enhanced and [Al/H] is underabundant by a 
factor of a few for the more plausible range of metallicities, — 1.4 < [Z/Zq] < —0.3. This pattern 



is similar to that seen in Galactic Halo stars of comparable metallicities ( Lauroesch et al. 1996 : 
Savage fc Sembach 1996|) . 



For Cloud B, we have measured iV(Aln), iV(Aliii), A(Cn), iV(Siiv) and an upper limit on 
A(Siii). As with Cloud A, the Lyman limit break provides the constraint of log A r (Hi) < 16.5 cm" 2 . 
Al is unusual in that its high-temperature dielectronic recombination rates are high and the ratio 
7V(A1ii)/A(A1iii) can increase with increasing U for log 2V(Hl) > 17 cm~ 2 flPetitjean et al. 1993 ). 
For smaller iV(Hi) this is not the case; a decreasing U yields an increasing A(Aln) as naively 
expected. Nonetheless, because of the unusual nature of reactions involving Al, we approached the 
modeling of Cloud B with caution. 

First, we considered the constraints that can be placed on the model cloud without 
considering the Al column densities (i.e. using only the Cll and Siiv detections and the Sin 
limit, we allowed Z, n#, and [Si/H] to vary). We found models consistent with the data for 
15.5 < logJV(Hi) < 16.75 cm" 2 . Smaller JV(Hi) values were unable to produce enough Siiv 
relative to Cn. For larger iV(Hl) values, U must be larger, corresponding to a smaller njj, and the 
model cloud was unrealistically large and Jeans unstable. For the acceptable range of iV(Hi), we 
found -0.3 < [Si/H] < -0.5. The metallicity ranged from [Z/Z ] = +0.2 for A(Hi) = 15.5 cm" 2 
to [Z/Zq] = —1.5 for A(Hl) = 16.75 cm~ 2 . Next, we considered what model clouds were consistent 
with the constraints provided by only the Alii and Aim column densities while allowing Z and 
uh to vary. Again, a consistent model cloud could be achieved for 15.0 < log iV(Hi) < 16.75 cm -2 . 
Systematically, the required metallicities were half to one dex larger than those required by the C 
and Si constraints alone. Naively, these two experiments imply [Al/Si]> 0, consistent with the 
striking strength of the Al profiles for Cloud B. 

When we ran CLOUDY in optimize mode, using all Si, C, and Al constraints and allowing 
Z, njj, [Si/H], and [Al/H] to vary, we found [Si/H] < and [Al/H] > were required for all 
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possible iV(Hi). Only the narrow range 15.9 < logiV(Hi) < 16.25 cm -2 yielded acceptable cloud 
models. For smaller values of iV(Hi), the ratios iV(Siiv)/iV(Siii) and iV(Aliii)/iV(Alii) could 
not be made consistent with the data (for the same U). The model results are illustrated in the 
lower panel of Figure [| Over the permitted range of iV(Hi), we have —0.55 < [Z/Z Q ] < —0.32, 
-2.5 <\ogU< -2, [Al/H] ~ +0.3, and -0.6 < [Si/H] < -0.2. 

Given the unusual conclusion we have reached that the Al in Cloud B is enhanced by a 
factor of two relative to solar, we now consider the validity of our assumption of photoionization 
equilibrium. Trapero et al. (1996) found the ratio of i\r(Aliii)/iV(Alii) out of equilibrium in 
the interstellar cloud toward the Bl V star 23 Ori. In that case, excited carbon implied a high 
density so that collisional ionization was applicable; the ratio of A r (Alin)/A r (Aln) implied a higher 
temperature than that measured from Doppler parameters, indicating that the gas had cooled 
faster than it could recombine. We obtained an unrestrictive upper limit on the density in Cloud B 
from an upper limit on iV(C*n). However, we have no reason to believe Cloud B suffers from such 
a non-equilibrium situation. Cloud B is consistent with photoionization in that the CLOUDY 
models are found to have kinetic temperatures consistent with those inferred from the Doppler 
parameters under the assumption of thermal broadening [such was not the case for the models of 
Trapero et al. (1996)]. 

Since the clouds are likely to be in a galaxy, it is possible that stellar sources could affect the 
shape of the ionizing radiation spectrum. However, the following basic argument applies. Both 
All and Sii have similar ionization potentials (6.0 eV and 8.2 eV, respectively) as do Alii and Sin 
(18.8 eV and 16.3 eV, respectively). Unless the spectrum of ionizing radiation has a very abrupt 
feature in these energy ranges, Alii and Sin will be similarly affected by a change in spectral 



shape. We ran test models and confirmed this with a star-forming galaxy spectrum (Bruzual & 
Chariot 1993|). 



4. Summary and Speculations 

Two kinematically adjacent, absorbing clouds in the same galaxy at v ~ km s _1 and 
~ 17 km s _1 have very different abundance patterns, if they are in photoionization equilibrium. 
Based upon CLOUDY models, Cloud A is inferred to have Si enhanced relative to a solar abundance 
pattern, and Al underabundant relative to solar. This abundance pattern is characteristic of Milky 
Way Halo stars flLauroesch et al. 1996| ; Savage fc Sembach 1996j ) which could have formed from 



gas like that observed in this z = 1.94 galaxy. We would expect such a pattern for many clouds in 
high redshift galaxies. Independent of the details, Cloud B has Al several times enhanced relative 
to Si (compared to the solar abundance pattern). This Al enhancement is highly unusual. If 
depletion onto dust grains is important we would expect a smaller Al/Si since Al depletes more 
readily than Si. To date, such an Al enhanced pattern has not been seen in absorption profiles 
looking through Milky Way interstellar gas ( Ravage fc Sembach 1996 ; [Snow et al. 1996] ; Spitzer & 



Fitzpatrick 1995; Fitzpatrick fc Spitzer 1994] ; |5pitzer fe Fitzpatrick 1995| ). 
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How could an enhancement of Al originate in a z = 1.94 cloud? In an attempt to find clues 
we note three astrophysical environments in which Al enhancement is observed: 1) in the stellar 
photospheres of the most metal poor globular clusters (Shetrone 199C; |Smith fe Kraft 1996[ ); 2) in 



the broad line regions of some AGNsQ ( [Shields 199% and 3) in the photospheres of Milky Way 



Bulge stars ( |McWilliam fc Rich 1991 ). Tnc 

common theme for the enhancement of Al in these 



three seemingly different environments is novae. The novae could produce the Al either directly in 
their ejecta ( jBmith fc Kraft 1996| ) or indirectly by providing magnesium isotopes to the ISM that 



later deposit onto stellar photospheres ( Langer fc Hoffman 1995 ). The 26 Mg and 27 Mg isotopes 



would then be converted to Al through proton capture in deep CNO convective mixing layers in 
metal poor stars ( Langer fc Hoffman 1995| ). More generally, does this suggest that a concentration 



of novae contributed to enhancing the Al in this cloud? In the globular cluster environment one 
key is to have a large number of stars that were formed coevally. Another key, which may apply 
to all three environments, is to have a potential well large enough to retain the gas. These may 
be prerequisites in our case also. Another possibility for excess Al production is a particular class 
of supernova for progenitors over a narrow mass range. The predicted amount of enhancement 



relative to other elements depends on the specific supernovae model adopted (|Nomoto et al. 1997[) , 
however only a small subclass of models would give Al/Si in agreement with that inferred for this 
cloud. 

Just how unusual is this cloud with a large Al/Si ratio? Such a pattern has not been seen 
in absorption along dozens of lines of sight toward Milky Way disk and halo stars. Most high 
redshift clouds do not have large Aim and Alii (relative to Sin). However, in the z = 2.14 



damped Lya absorber in Q0528 — 251 (Lu et al. 1996|), we have identified a single outlying cloud 



(at ~ —100 km s~ x ) for which the equivalent width of the Alii A1671 transition is larger than 
that of the Sin A1527 transition. The other clouds in the same system clearly have the opposite 
ratio. Finding more examples and establishing similarities between Al-rich environments are 
logical next steps in diagnosing the origin of this abundance pattern at high redshift and perhaps 
understanding the anomalous enhancements seen in metal-poor globular cluster and Bulge stars. 
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4 We note that the absorber lies 30,000 km s 1 from the emission redshift of the QSO, which is not a BAL AGN. 
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Table 1. Column Densities of Key Transitions 



Transition 


Cloud A 


Cloud B 


log(N cm 2 ) 


6 (km s 1 ) 


log(N cm 2 ) 


6 (km s 1 ) 


Alii 1670.787 


< 11.13 a 




11.85 ±0.04 


4.91 ±0.60 


Aim 1862.790 


< 11.79 a 




12.06 ±0.10 


5.17 ± 1.71 


Cn 1334.532 


13.11 ±0.11 


4.80 ±2.02 


13.45 ± 0.07 


6.91 ± 1.44 


Sin 1526.707 


12.42 ±0.11 


1.84 ±1.88 


< 12.29 b 




Siiv 1393.755 


12.72 ± 0.04 


4.70 ± 0.74 


13.11 ±0.03 


6.50 ±0.53 


Si iv 1402.770 


12.72 ± 0.04 


4.70 ± 0.74 


13.11 ±0.03 


6.50 ±0.53 



a Upper limit for assumed 6 = 5 km s 1 , based upon the Cll A1334 transition 
b Upper limit for assumed 6 = 5 km s _1 , based upon the Alii A1670 transition 
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Fig. 1. — The HIRES /Keck absorption profiles of the z = 1.94 Civ system in the spectrum of 
Q1222+228. The two clouds studied here, A and B, are marked by the vertical ticks at v ~ km s _1 
and v ~ 17 km s _1 , respectively. The cross mark at v = 115 km s _1 near the Cn A1335 transition 
flags an unidentified absorption feature. 
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Fig. 2. — The [Al/H] and [Si/H] abundance ratios (relative to solar) as a function of log iV(Hi) for 
clouds A (upper panel) and B (lower panel). The corresponding metallicities, [Z/Zq], are given by 
the upper horizontal axes. The unshaded regions give the range of models that are consistent with 
the observed data (see text). The long-dash curves are the ionization parameters, logU, which are 
given by the right hand vertical axes. The iV(Hl) at which the model cloud sizes exceed 10 kpc 
and become Jeans unstable are marked by the vertical short-dash lines. 



